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<^ ■ Abstract 

, We investigate the shell structure of bubble nuclei in simple phenomenological 

shell models and study their binding energy as a function of the radii and of the 
number of neutron and protons using Strutinsky's method. Shell effects come about, 
on the one hand, by the high degeneracy of levels with large angular momentum 
\^ \ and, on the other, by the big energy gaps between states with a different number of 

■ radial nodes. Shell energies down to -40 MeV are shown to occur for certain magic 

. nuclei. Estimates demonstrate that the calculated shell effects for certain magic 

I numbers of constituents are probably large enough to produce stability against fis- 

0^ ' sion, a-, and /?-decay. No bubble solutions are found for mass number A < 450. 

PACS numberes: 21.10.Sf,24.10.Nz,47.20.Dr,47.55.DZ 

F; ; 1 Introduction 

> . 

^ . The possibility that nuclei could exist in the form of a (spherical) bubble or in the form of 
a toroid has been pointed out long ago [1,2]. Within the liquid drop model (LDM) nuclei 
of these shapes turn out to be unstable with respect to deformations. Shell-effects may, 
however, stabilize such nuclei against deformation. 

In a series of papers [3], C.Y. Wong investigated shell-effects for toroidal and bubble- 
shaped nuclei using Strutinsky's shell correction method. He restricted his attention to 
known nuclei near the valley of /3-stability and found that for certain doubly magic nuclei 
(|Q°if(y'i2o, Ceso) Spherical bubble solutions with a very small inner radius (ratio of 
inner to outer radius ~0.07) turned out to be the ground state. Indications that bubble 
solutions might exist were also found in mean field calculations [4] and for stellar matter 
at finite temperature [5]. More recently, Moretto et al. [6] showed in a classical model 
that LD-bubbles at finite temperature may be stabilized by an internal vapor pressure. 

In the present paper, we study shell effects for nuclear bubbles in a broad range of 
neutron (A^) and proton (Z) numbers extending considerably beyond the known nuclei. 
As C.Y. Wong we make use of Strutinsky's shell correction method [7]. We show that the 

*We dedicate this paper to Prof. Richard Lemmer on the occasion of his 65*'* birthday. 
tOn leave on absence from University M.C.S. in Lublin 
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shell energy may become as large as -40 MeV for certain magic numbers of the nuclear 
constituents and that nuclear bubbles may thus become stable or very long-lived against 
fission and other decay modes. 



2 Theory 

In Strutinsky's method [7], the total binding energy E of a. nucleus (neutron number N, 
proton number Z, nucleon number A) of given shape is given as a sum of the liquid drop 
(LD) energy £^ld(-^, Z) and the shell correction energy (5£'sheii 

E{N, Z) = ELD(iV, Z) + 5E,hell(iV, Z) . (1) 

The shell correction energy has the well-known form [7] 

SE^heii = ^(iu Sn^ , (2) 

where Ci, are the single particle energies and 5ni, is the difference between the occupation 
numbers rii, — 9o{eF — Si) in the shell-model ground state and a smooth distribution 

Sn^ = n^-ny. (3) 

The smooth occupation pattern is defined in the usual way as a functional of a smooth 
level distribution [7]. 

We use the Strutinsky method to study the total energy and especially the shell cor- 
rection energy of spherical nuclear bubbles. The single particle energies e^, as well as the 
LD energy depend on the inner {R2) and outer {Ri) radius of the bubble nucleus. Adopt- 
ing the conventional saturation condition that the volume of the LD remains constant 
independently of its shape, the two radii are related by the condition 

where Rq = tqA^/^ is the radius of a compact spherical nucleus of the same mass. We 
describe the shape of the bubble nucleus either by the dimensionless radii 

^^1(2) (5) 
or by the ratio / between the volume of the hole and the volume of the entire bubble 



R. 



The difference AE'ld between the energy of the LD-bubble and the energy of the corre- 
sponding compact spherical LD is given by 

AE'ld = Es{Ri, R2) + Ech{Ri, R2) — Es{Ro) — Ech{Ro) , (7) 
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where Eg and Ec\, are the surface and Coulomb energies of the bubble and the com- 
pact configuration, rcsp. Measuring Ai^LD in units of the surface energy Es{Rq) and 
substituting the explicit form of the different energies, we find 



E := 



LD 



Es{Rq) 



^vl{v2)+vl-l + 2Xo 



VI{V2) + - ^Vlvl{v2) 



where 



vi{v2) := [l + vlY^ 
and Xq is the conventional fissility parameter 



-E^cb(-Ro) 



Z^/A 



2Es{Ro) {ZyAUt 



with 



3eg 



(8) 



(9) 



(10) 



(11) 



The quantities a, tq, and Cq are the surface constant, the radius parameter, and the 
elementary charge, resp.. Beside the trivial solution — Q (compact spherical nucleus) 
the condition of stationarity 

1^ = (12) 

has real solutions ^2(^0) only for Xq > 2.02. They represent a multivalued function which 
is shown in the l.h.s. of Fig. 1 for a limited range of values V2{Xq). For given Xq > 2.02 
the smallest solution V2{Xq) < 0.4 corresponds to the maximum of the potential barrier, 
whereas the next solution in magnitude V2{Xq) > 0.4 corresponds to a minimum of the 
energy S for a bubble of the reduced radius V2{Xq). 

On the r.h.s. of Fig. 1 we show the energy change S by bubble formation as a function 
of V2. For Xq > 2.2, the bubble solution is seen to correspond to a lower energy than 
the one of a compact sphere (^2 = 0). The barrier between the bubble solution and the 
compact spherical one occurs for reduced inner radii in the interval < V2 < OA. In units 
of Es{Ro) the barrier heights are seen to be at most at S = 0.025. The unit Eg = AnaR'^ 
is, however, pretty large. For A = 580, Z = 280 and reasonable values of the parameters 
To and a, we have Es = 1291 MeV and Xq = 2.62. The corresponding energies of the 
barrier and the minimum are ~ 20 MeV and ~ —140 MeV. It is not meaningful to pursue 
the solutions V2 of Eq. (12) to values of V2 and Xq much larger than the ones given in Fig. 
1, because for increasing Xq the charge density of the bubble turns out to be unphysically 
large and the diameter d = Ri — R2 oi the bubble layer becomes unphysically small. 

We emphasize that the bubble solutions obtained in the LDM are not stable with 
respect to deformations [3] in the same way as the compact spherical liquid drops are 
not stable agains fission for fissibility parameters Xq > 1. Nevertheless, stabihty can in 
principle be produced by shell effects. 
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As two extreme and simple cases of nuclear single particle potentials we consider the 
shell effects in an infinite square well and in a harmonic oscillator: 



y(r) = | for i?,<r<i?i 

^ ' y +00 otherwise ^ ' 

V{T)^-V, + ^{r-Rf. (14) 

The depth Vq > has no influence on the shell correction energy and can thus be put 
equal to zero. The center of the oscillator potential is chosen to be 

R=?^±^. (15) 

The oscillator frequency u is chosen in such a way that the RMS deviation from the sphere 
of radius R is the same when calculated with the shell model wave function and with the 
LD-density 

((r - ^)')sM = ((r - ^))ld . (16) 

We have to add a spin-orbit term to the central potential (13) or (14). We use the 
conventional form in the Skyrme approach (see Eq. (5.103) in Ref. [8]) 

Vso = Vso(r)r-l', (17) 

p(r) is the nuclear density distribution. We choose the value Wq = 120 MeV fm^ given in 
Ref. [8] and neglect the isospin dependence which is recently under debate [9]. 

For the total density p(r) which appears in the expression for the spin-orbit potential 
Vso('") (Eq- (18)) we used the following ansatz: 



( r-R 

1-2 



Ri — R2 



(19) 



where po=0.17 fm~^ and R is given by Eq. (15). For simplicity we have assumed that the 
proton and neutron densities Pp and Pn are proportional to Z or N respectively. 

It is seen from (18) that the sign of the spin-orbit potential is opposite in the inner 
and outer surface region of the bubble nucleus. Consequently, the magnitude of the spin- 
orbit splitting is smaller for bubble nuclei than for normal ones. We, therefore, treat the 
spin-orbit term in perturbation theory. The eigenenergies Cnij of the s.p. Hamiltonian 
including the spin-orbit potential (17) and the unperturbed eigenvalues Sni are related by 
the equation 

enlj = ^nl + i'ipnljmlVsol'ipnljm) ' + 1) " + ~ i) ' ^^^^ 
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where (n — 1) represents the number of radial nodes (not counting zeros at r = 0, oo) and 
l,j,m are the orbital and total angular momentum and its projection, resp.. The mean 
value {ipnijmlVsoli'nijm) dcpcuds only on the unperturbed s.p. density 

Pnlir) = i^tl.^i^i^nlU^ = ^ . (21) 

The functions Uni{r) satisfy the radial Schrodinger equation with eigenvalue Sni- For 
the infinite square well the eigenvalues and eigenfunctions are obtained by numerically 
satisfying the boundary conditions and for the harmonic oscillator (14) we used the WKB 
approximation. 



3 Results and Discussion 

We complement the results obtained in the pure LDM (see Fig. 1) by Fig. 2 which shows 
lines of equal LD-binding energy per particle {Eu:,{N, ^ function of and Z. 

We note that the volume term of the LDM is present in this figure since we do not subtract 
the energy of the compact spherical LD. The LD parameters are taken from Ref. |Ty . 



Each point on the equipotential lines corresponds to a spherical bubble solution as given 
by the LDM, which is unstable with respect to deformations. It should be noted that the 
largest energy gains per particle occur for nucleon numbers 1200 < A < 2000 and the 
corresponding proton numbers 325 < Z < 400. Some of the isobaric lines are cut twice 
by the same equi-energy line. If the binding energy per particle in between these 2 points 
is smaller, a /3-stable isobar lies somewhere on this section. Of course, this consideration 
is thwarted by the fact, that the LD-bubbles are all unstable against fission. In the LDM, 
the (spherical) bubble solutions are saddle points, not minima. 

The LD results (and consequently also the result on the total binding energy) depend 
sensitively on the value of the surface constant a. This surface constant contains a poorly 
known isospin-dependent part. In our calculations we adopted the isospin dependence 



given by Myers and Swiatecki [|T0[ who write the surface energy of a spherical LD in the 
form 

'N-Z' 



Es{Ro) = 17.9439 MeV • A^^'^ 



1 - 1.7826 



A 



[22) 



In Fig. 3 we show the spectrum of single particle levels for the shifted infinite square 
well with a spin-orbit term. The levels are shown as a function of the hole fraction 
parameter / (see Eq. (6)). As / increases, the diameter d = R2 — Ri of the bubble layer 
decreases. Increasing the number of radial nodes for given orbital angular momentum I 
thus costs an energy which rises steeply as a function of /. Augmenting a given /-value by 
1 for given n and given parameter / implies an increase of the centrifugal energy, which 
is the larger, the larger the /-value. Magic numbers come about by the interplay between 
levels Cnij with n > 2, which rise rapidly as a function of / and levels euj, which depend 
more gently on / with a tendency to decrease with / due to the diminishing centrifugal 
energy. 
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In Fig. 4, we display a corresponding level spectrum for the harmonic oscillator 
potential. The general trends are the same as for the infinite square well, but the magic 
numbers for the same values / differ for the square well and oscillator. This is not 
surprising. A Saxon- Woods form centered around R would be more realistic and the 
magic numbers for this choice would be expected to lie between the limits given by the 
infinite square well and the oscillator. 

For both the potentials we observe that the spin-orbit splitting is often reversed as 
compared to the case of normal nuclei. The reason was already given in Section 2. The 
shell correction energy 5i?sheii is shown as a function of N (or Z) at a given value of 
/ = 0.28 for the infinite square well in Fig. 5 and for the harmonic oscillator in Fig. 6. 
Please note the different energy scales in case of the square well and the oscillator. The 
eigenenergies of the square well scale with A^"^^^ which is taken into account by the energy 
unit. No such simple scaling property exists for the oscillator. It is seen from the Figs. 
5 and 6 that the shell energy may produce energy gains up to -20 MeV for one sort of 
particles. Of course, double magic shell closures can only occur, if the two magic numbers 
correspond to the same /-values. 

In Fig. 7 we show lines of constant shell correction energy 5-Eshcii in the (A^, Z)-plane. 
Each point corresponds to a minimum of the total energy as a function of /. One should 
notice that for this case of relatively light bubbles, the shell effects are especially large for 
almost symmetrical nuclear composition. 

In Fig. 8 we display lines of equal energy gain by formation of a bubble. As a reference 
we use the energy of a spherical LD of the same numbers of neutrons and protons. It is 
seen that the gain in binding energy may amount to several 100 MeV. 

We have still to deal with the crucial question of stability of the bubble solutions with 
respect to shape deformations. Calculating Strutinsky's shell correction energy (2) for a 
deformed bubble nucleus implies that we find the eigenvalues of a Schrodinger equation in 3 
dimensions or at least (for axial symmetry) in 2 dimensions. Without solving this technical 
problem we may obtain a valuable insight by describing the deformation dependence of 



the shell correction energy as suggested by Myers and Swiatecki [11]. In Fig. 9 we display 



the dramatic effect of the shell correction on the binding energy of a bubble nucleus as 
a function of the quadrupole deformation. The total LD binding energy of the spherical 
bubble nucleus is put equal to zero. The LD-part of the energy decreases monotonically 
as 1/32 1 increases. Adding the shell energy (5i?sheii with Swiatecki's /?2- dependence produces 
a valley of about -30 MeV depth. The barriers on both sides of the minimum are about at 
an energy of -3.9 MeV which leaves us with a barrier height of about 25 MeV. This order 
of magnitude of the fission barrier implies an almost vanishing probability for spontaneous 
fission. 

What about the other decay modes? /?-decay will imply that the system moves along 
isobaric lines in the (A^, Z)-plot towards lower energy. Fig. 8 is particularly instructive in 
this respect: If a given isobaric straight line cuts a line of constant binding energy twice 
and if, at the same time, the binding energy between the two points of intersection is 
lower, a /3-stable nuclear bubble lies in between. In Fig. 8 this happens to be the case 
for A > 600. If, on the other hand, the energy keeps on lowering along a given isobaric 
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direction, the /3-decays may finally lead to a nuclear composition where bubbles cease to 
exist. Thus there are cases where bubble nuclei may disintegrate by a series of /3-decays 
and others where the /9-decays make them approach a stable composition. 

The a-dccay, which limits the lifetime of the presently known superheavy nuclei, cer- 
tainly may also limit the lifetime of bubble nuclei. The penetrability of the Coulomb 
barrier for an a-particle depends exponentially on the Coulomb potential at r = Ri, 
which has the vahic 2Zeo//?i. The higher the Coulomb potential, the lower the a-decay 
probability. For A = 700, Z = 270, and i?i = 11 fm one finds 2Zel/Ri ^ 70 MeV. We 
estimate the Q- value to be 



Thus the a-particle has to tunnel through a Coulomb barrier whose top is ~50 MeV 
above the tunneling energy. Empirically, the a-decay of heavy actinide nuclei is found to 
be vanishingly small whenever the emitted a-particle has an energy Qa < 4 MeV, i.e. if 
the energy difference 



This estimate implies that the a-decay probability of typical bubble nuclei is very small 
due to the especially large Coulomb barrier. Of course, much more careful investigations 
of the decay modes have to be performed. We may, however, safely say that certain 
bubble nuclei are expected to have a practically infinite lifetime. 

It is intriguing to imagine that stable or at least very long-lived bubble nuclei may 
exist. Their properties would be a fascinating subject of research. Unfortunately, the 
masses and charges of the best candidates for bubble structure are so high that there is 
no hope to produce them. More careful work with the full Hartree-Bogoliubov theory is of 
course necessary, especially for determining the lower limits of mass and charge numbers 
of these objects. Even if bubble nuclei can never be made in an earthly laboratory, they 
might play a role in neutron stars [5]. 

Finally, there may be bubble structures for mesoscopic systems consisting of some 
1000 atoms. This was already conjectured in Ref. [6]. We beheve that it may be also 
of interest to investigate whether shell effects favour the bubble topology for certain 
mesoscopic systems. 

Acknowledgments 

K. Pomorski is grateful to the DFG for granting a guest professor position which 
greatly facihtated the present research. Prior to the 1^* April 1997, this research was 
also modestly supported by the BMBF which is gratefully acknowledged. This work was 
also partially financed by the Polish Committee of Scientific Research under contract No. 
2P03B 049 09. K.D. is grateful for numerous stimulating e-mail messages by M. Weiss. 



g« := E,,,{Z, N) - [E,,,{Z - 2, AT - 2) - 28 MeV] 20MeV . 




7 



References 

[I] P.J. Siemens and H.A. Bethe, Phys. Rev. Lett. 18, 704 (1967). 
[2] J. A. Wheeler, Notebook, unpublished . 

[3] C.Y. Wong, Ann. Phys. (N.Y.) 77, 279 (1973) (with reference to earlier work). 

[4] P. Bonche, P. H. Heenen, M. Weiss, private communication. 

[5] P. Bonche and D. Vautherin, Nucl. Phys. A372, 496 (1981). 

[6] L.G. Moretto, K. Tso and G.J. Wozniak, Phys. Rev. Lett. 78, 824 (1997). 

[7] M. Brack, J. Damgaard, A.S. Jensen, H.C. Pauli, V.M. Strutinsky and C.Y. Wong, 
Rev. Mod. Phys. 44, 320 (1972). 

[8] P.Ring and P. Schuck, "The Nuclear Many Body Problem", New York 1980, Springer 
Verlag. 

[9] P.G. Reinhard and H. Flocard, Nucl. Phys. A584, 467 (1995) and ref. therein. 
[10] W.D. Myers and W.J. Swiatecki, Ark. Fys. 36, 343 (1966). 

[II] W.D. Myers and W.J. Swiatecki, Nucl. Phys. 81, 1 (1966). 



8 



Figures captions. 



1. The condition (12) defines a function V2{Xq) wliicli is sliown in tlie l.li.s. of the 
figure. The r.h.s. shows the difference £ = AEi^d/ Es{Ro) of the LD energy of a 
bubble of reduced inner radius f 2 and the LD energy of a compact spherical nucleus 
of radius Rq. 

2. Lines of constant LD-binding energy per particle as a function of {N,Z). Each 
point on the lines corresponds to a bubble solution within the LDM. The straight 
lines in the plot represent isobars. 

3. Level scheme as a function of the hole fraction / = {R2/ RiY for the infinite square 
well plus the spin orbit term (Eqn. 17-19). For 3 values of / (0.12, 0.24, 0.28) 
numbers are written just above certain s.p. energies whenever the distance {cv+i — 
e^) exceeds 1.5 energy units. Following the conventional spectroscopic notation the 
1- values of the levels are given by the letters of the alphabet. The energy unit used 
takes into account that the eigenvalues of the infinite square well scale with A"^/^. 

4. As Fig. 5, but for the harmonic oscillator potential. The oscillator frequency uj was 
determined as a function of / for y4=500 using relation (16). 

5. Shell correction energy as a function of Z (or A^) for /=0.28 for the infinite square 
well. 

6. As in Fig. 7, but for the harmonic oscillator potential. The frequency of the 
oscillator was determined by the relation (16) for the mass number ^4=500. 

7. Lines of constant shell energy 6Esheii in the (A^, Z) plane. SEsheii was calculated for 
the infinite square well. Each point on the equi-energy lines is evaluated for the hole 
fraction parameter / corresponding to the minimum of the total energy. 

8. Lines of constant energy gain (Ai^LD + SEsheii) with respect to the energy of a 
compact spherical LD. The shell correction energy was calculated with the infinite 
square well potential. 

9. LD-energy (dashed line), shell correction energy (dotted line) and total energy (solid 
line) as a function of the quadrupole deformation (32 of the outer bubble surface 5*1. 
The liquid drop energy was minimized with respect of the deformations jS^ and (3q 
of the outer surface and j32 to of the inner surface. The deformation dependence 



of the shell correction was taken from Ref. Ill 
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